On the integrability of halo dipoles in gravity 
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Abstract 



We stress that halo dipole components are nontrivial in core-halo systems 

in both Newton's gravity and General Relativity. To this end, we extend a 

recent exact relativistic model to include also a halo dipole component. Next, 

we consider orbits evolving in the inner vacuum between a monopolar core and 

a pure halo dipole and find that, while the Newtonian dynamics is integrable, 

its relativistic counterpart is chaotic. This shows that chaoticity due only to 

halo dipoles is an intrinsic relativistic gravitational effect. 
PACS numbers: 04.20.Jb, 05.45.+b, 95.10.Fh, 95.30.Sf 



Typeset using REVT^ 



mail: vieira@ime.unicamp.br 
mail: letelier@ime.unicamp.br 



1 



We report in this Letter a simple situation where the contact between Newton's gravity 
and General Relativity (GR) allows for a clear distinction between their predictions about 
chaotic behavior of orbit dynamics. 

We considered in [|I| how does GR deal with a monopolar core (at first approximation), 
which may be a neutron star, a galactic nucleus or even a black hole, surrounded by a distant, 
massive halo of dust. We exhibited an exact, static, axially symmetric model with a halo 
made of quadrupoles and octopoles of arbitrary strengths. In constructing this solution, we 
were also guided by its formal resemblance to the classical Henon-Heiles potential in the 
Newtonian limit. 

Here, we are mainly interested in the dynamical rather than formal comparison between 
Newtonian core-halo systems and its relativistic counterpart. Firstly, we state the problem 
according to Newton's theory: let the coordinate origin stay at the mass center of the 
monopolar core, z be the symmetry axis of the core-halo system and D be the region 
between the minimum and the maximum spheres centered at the origin that isolate the 
inner vacuum from the core and the halo. We have to solve Laplace's equation in D for the 
axially symmetric Newtonian potential. By using standard spherical multipolar expansion, 
we arrive at the following gravitational potential felt by test particles evolving in D: 

$,v = $o---2^2 + f (22'-r') + ^(2^'-3zr2) + ... . (1) 

$0 is a constant (put here only for the sake of completeness) and = r"^ + z"^ = + + z^ 
where x,y,z are the usual Cartesian coordinates. $Ar and all quantities in it are written 
in convenient nondimensional units. V, Q and O are respectively the dipole, quadrupole 
and octopole halo strengths given by the appropriate normalized integrals over the halo 
(as a source) in the Newtonian approach. It is worth to stress that, quite apart from the 
nomenclature, the halo multipoles have opposite behavior in D from the usual decreasing- 
with-distance core multipoles. In particular, by merely translating the origin to the core 
mass center we can (and in fact we do) eliminate the usual core dipole in D, so that the next 
higher core contribution should be the usual quadrupolar one (which will not be considered 
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here). On the other hand, nothing so simple can ehminate the halo dipole component. As 
it is easily seen, the permanency of the later is related to the halo mass distribution with 
respect to the plane z = 0. The only way of eliminating halo dipoles from (|1]) should be 
passing to a noninertial reference frame with uniform acceleration V in the z-direction. 
In this case, we should have to properly treat this problem as a relativistic one, since it 
extrapolates the Galilean gauge invariance of the Newtonian physics. 

Now, we consider bounded orbits under the Newtonian potential If we switch the 
halo quadrupole and octopole off (Q = O = 0), thus the remaining potential describing 
the core plus the pure halo dipole is integrable, in particular being separable in parabolic 
coordinates 0. The hamiltonian H as well as the additional constant of motion C in 
involution with the hamiltonian are 

C = rprP, - zpl - + - - ^r^, (2) 

where ^ is the nondimensional conserved angular momentum associated to the axial sym- 
metry. Only to illustrate the integrable dynamics, we present in Fig. 1(a) a typical z = Q 
Poincare section for this case. Contrasting with this, we mention that either the Newtonian 
halo quadrupole (Q 7^ 0) or octopole {O 7^ 0) does produce chaos, which will be shown in 
an extended forthcoming paper. 

Since the Newtonian halo dipole is nontrivial in appearing quite naturally in (P, 
we ask for its relativistic counterpart. In fact, we succeeded in extending our previous 
relativistic solution to add a dipole term to the halo. Starting from Weyl's metric and 
following the formulation presented therein, the solution depends upon two functions v{u, v) 
and 7(m, f ), the last being obtained by quadrature from the former one. Here, it is sufficient 
to write down only the main function v{u^v) (see for notation and details): 

v{u, v) = aMu) - VPi{u)Pi{v) + {2/'i)QP2{u)P2{v) + {2/b)OP:,{u)P^{v). (3) 

This equation should be compared with eq. (3) of that Reference. Essentially, the first term 

3 



with Oo = —1 describes the monopolar core (that can be switched off by putting oq = 0), the 
remaining ones being the corresponding muhipoles originated from the halo. Qq and P„ are 
the corresponding Legendre's functions in standard notation and u, v are prolate spheroidal 
coordinates constructed from Weyl's ones. 

What is the Newtonian limit of (^? The basic steps to obtain it are: i) The relativistic 
solution is valid in the full vacuum between the core and the halo, ii) We assume that there 
exist a region D in the intermediate vacuum where the conditions of weak gravitational 
field and slow motion of test particles occur. Then, Eintein's equations reduce in D to 
Laplace's equation for the Newtonian potential $Ar, which relates to the metric through the 
temporal component gu = I + (2/c^)$7v (the remaining spatial components of the metric 
being irrelevant to this approximation), iii) The Schwarzschild coordinates at which we 
finally come back with the solution approximate to the usual Euclidean spherical ones plus 
time in D. We consistently assume that this region is far from the core's event horizon. Now, 
by applying all this to the relativistic solution (^ and expanding gtt to the first order in P, 
Q and O we arrive, after the conversion to nondimensional units, exactly at the potential 
$iv given by ([I|) (apart from the irrelevant constant $o)- The implicit assumption here is 
that D and the Newtonian region D above have a nonempty intersection at least. We would 
like to stress the nontrivial character of halo dipoles also in the relativistic context, as they 
contribute to the Kretschmann's scalar R^^^"'^ R^p^s even if oq = Q = C = 0. Moreover, they 
are smooth, however exact, relativistic counterparts of Newtonian halo dipoles, in the sense 
that their first contribution to the Riemann tensor is of second order in V. 

Next, we compare bounded geodesies generated by (|^) for a pure halo dipole (Q = C = 0) 
with the preceding Newtonian case. A typical Poincare section of relativistic orbits is shown 
in Fig. 1(b) where vast chaotic layers are seen in the geodesic motion. This should be 
contrasted with the full integrability of the corresponding Newtonian figure 1(a). Since 
Poincare sections are global, gauge invariant tools to detect chaos in bounded systems, this 
proves that chaoticity in the vacuum between a monopolar core and a purely dipolar halo 
of dust is an intrinsic relativistic gravitational effect. 
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An extended numerical study of both newtonian and relativistic well as the 

details of the full relativistic solution (||) will be presented in an extended forthcoming 
paper. 

The authors thank CNPq and FAPESP for financial support. 
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FIGURES 

FIG. 1. Poincare sections of orbits in the vacuum between a monopolar core and a pure halo 
dipole (2 = = 0) through the plane z = (for pz > 0): (a) the Newtonian case in Euclidean 
cilyndrical coordinates for the potential (||) with nondimensional orbit parameters E = —0.03 
(energy), £ = 3.8 (angular momentum) and T> = 5. x 10^^ (halo dipole strength), and (b) the rela- 
tivistic case for timelike geodesic motion in Weyl's coordinates for the solution (^) with E = 0.975, 
i = 3.8 and V = 3. x 10~^. Note that the relativistic energy E = 0.975 corresponds in the 
Newtonian limit to the energy E = —0.025. 
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